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We study the dynamics of adatoms in a model of vicinal 
(11m) fee metal surfaces. We examine the role of different dif- 
fusion mechanisms and their implications to surface growth. 
In particular, we study the effect of steps and kinks on adatom 
dynamics. We show that the existence of kinks is crucially 
important for adatom motion along and across steps. Our 
results are in agreement with recent experiments on Cu(lOO) 
and Cu(l,l,19) surfaces. The results also suggest that for 
some metals exotic diffusion mechanisms may be important 
for mass transport across the steps. 

PACS numbers: 61.50.Cj, 68.35. Fx, 68.55.Bd 



Diffusion of adatoms on solid surfaces is an extensively 
studied subject In particular, adatom dynamics on 
vicinal metal and semiconductor surfaces has important 
implications to surface growth under non-equilibrium 
conditions [||J|]. However, barring a few special cases, 
the atomistic details of diffusion processes near steps and 
kinks are not known. Current experimental techniques 
are now able to yield atomistic information about adatom 
dynamics near steps and kinks [Q , and growth of surfaces 
. Clearly, careful microscopic calculations are needed 
to understand these phenomena. 

Our aim in this Letter is to study models of surfaces 
of fee metals vicinal to the (100) plane. The open struc- 
ture of the (100) facets can be expected to give rise to 
some unconventional diffusion processes that are not seen 
for instance on stepped surfaces with fcc(lll) terraces. 
First, we want to identify the various microscopic mecha- 
nisms relevant to self-diffusion near steps and kinks. Sec- 
ond, we shall discuss the implications of our results to the 
morphological stability of these surfaces under growth , 
and suggest a phenomenological model for step growth. 
Our results are consistent with experiments on Cu(lOO) 
§ and Cu(l,l,19) surfaces. 

The geometric structure of an fcc(119) surface is shown 
in Fig. 0. An ideal fcc(llTO) facet, with odd to > 1, 
consists of (100) terraces of width (m — l)r„„/2 sepa- 
rated by (111) steps of height r„„/\/2, where r„„ is the 
distance between nearest neighbor atoms. Due to the 
geometry, only one kind of steps (of monolayer height) 
with close-packed edges exist on these surfaces. The 
metallic interactions between atoms in our model are de- 
rived from the semi-empirical Effective Medium Theory 
(EMT). The formaUsm of EMT is presented in Ref. 0, 



and a description of the implementation for molecular 
dynamics (MD) simulations of the present work can be 
found in Ref. |8|. In the case of copper, EMT has been 
shown to give a reasonably accurate quantitative descrip- 
tion of many different surface phenomena j? - 1 1 1 , which 
motivates its use for the present case. 

We shall divide the discussion of the microscopic mech- 
anisms near step edges into three parts: standard hop- 
ping events (denoted by H), exchange and other exotic 
mechanisms (X) , and the effect of kinks on diffusion near 
and across step edges (K). For each mechanism M the 
activation barrier is denoted by i?M, and that of the re- 
versed process by . 

In Fig. ||(a) we show a contour plot of the adiabatic 
energy surface E{x, y) experienced by a single adatom on 
the Cu(119) surface at zero temperature. The potential 
across the terrace is shown in Fig. ||(b), indicating the 
activation energy for diffusion on the terrace £^a, and 
the height of the Schwoebel step barrier E'^^ |l^. It is 
evident that the barrier in the x direction is appreciably 
modulated only at the immediate vicinity of the steps 
p3[ . We have verified this for the Cu(l,l,15) surface 
also. 

Activation energies for simple hopping mechanisms on 
surfaces of several fee metals with different orientations, 
as given by EMT, have been extensively tabulated in 
Ref. (1^. Our results for copper are fully consistent. 
The barrier height for a single jump on a flat terrace 
far from step edges is found to be i?A — 0.399 eV, and 
that for diffusion of a vacancy in the first layer of the 
terrace is Ey = 0.473 eV. The lowest barrier is that of 
an adatom diffusing along a straight step edge, and has a 
value i?H2 = 0.258 eV. As expected, on Cu(IIto) surfaces 
we find i?H2 < Ea< ^^hT < -^Hi- The inequality i?H2 < 
i?A is consistent with experimental results ||l^ . 

In addition to ground-state calculations we have per- 
formed extensive MD simulations ||l^ to identify possi- 
ble exotic diffusion mechanisms Jl6| and to study entropic 
contributions to the rates [ pT| . A well-known mechanism 
for step crossing, the replacement of an edge atom by an 
adatom from the terrace (mechanism XI) is observed in 
our simulations. In our model the activation energies for 
hopping and the simple exchange across the step edge are 
approximately equal: i?Hi ~ £^xi |p8 19|. We have also 



found more complicated mechanisms for step crossing. In 
Fig. ^we show two examples: a "coherent" chain transfer 
and an atom-by-atom replacement mechanism (vacancy 
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diffusion). A possible explanation for the first mechanism 
is the local close-packed-like order of the surface atoms 
in the second configuration of Fig. ||(a) , which may lower 
the local free energy. This is obviously an effect char- 
acteristic to stepped surfaces with fcc(lOO) terraces: e.g. 
on an fcc(lll) terrace the atomic rows are more densely 
packed and cannot easily slide with respect to each other. 
The second process shown in Fig. |3|(b) can be described 
as a "popping up" of a surface atom onto the step edge 
and the diffusion of a vacancy towards the descending 
step. By reaching the step, the vacancy turns into a hole 
or a pair of kinks at the step edge, which can then be 
filled by a surface atom or an adatom from the terrace 
below. Repeating this procedure, e.g. under an exter- 
nal field driving the atoms into the negative x direction, 
can result in mass transfer across steps which can then 
enhance growth instability . 

The activation energy of the first stage of the pro- 
cess shown in Fig. ||(b), i.e. the pop-up of a surface 
atom, was reduced by the existence of a kink at the 
step edge. For the processes with and without a kink 
we find Ek3 < £^X3, respectively (cf. Fig. ||). Indeed, 
the effect of kinks on the energetics of diffusion processes 
near step edges seems to be crucial. In Fig. || we show 
the most important diffusion routes near straight and 
kinked step edges. It turns out that the hopping of a 
single adatom across the step edge in the vicinity of a 
kink site (Kl) is not much more favorable than climb- 
ing across a straight edge (HI). On the other hand, 
activation barrier for the escape of an adatom from a 
kink (K4) is higher than Ea, while going "around the 
corner" along a kinked step edge (K2) is even more ex- 
pensive, i.e. £'h2 < Ea < Eka < £^k2- From exper- 
iments on Cu(l,l,19), the activation energy for mass 
transport along kinked step edges was determined to be 
10300 ± 1630 K 0], which was assumed to be due to 
the process K4. For K4 we obtain 6011 K in agree- 
ment with Ref. |lj]. However, if we assume that K2 
is the rate-limiting process we get an activation energy 
of £'k2 — 9075 K which is within the experimental error 
bars. 

The activation barriers for the processes shown in 
Fig. ^ for Cu(llr7i) are summarized in Table |. For 
any mechanism of step crossing in either direction, the 
barrier height is found to be well above Ea, i.e. a clear 
Schwoebel barrier exists. Our results thus indicate that 
under growth conditions the currents should go upward 
which makes the (100) surface unstable This is con- 
sistent with the experimental results of Ref. on growth 
of the Cu(lOO) surface. 

In the case of our copper model, the activation barriers 
for HI, XI, and K3 across the step edge are almost the 
same. This means that a finite density of kinks promotes 
step crossing. Due to small differences in activation en- 
ergy, the relative occurrence of the different mechanisms 
at low temperatures is expected to be strongly model 



and material dependent. In particular, our results sug- 
gest the possibility that for (11m) surfaces of other fee 
metals, exotic mechanisms such as K3 could play a more 
important role in mass transport across step edges. In 
addition, some processes such as X2 are influenced by 
the step density. More systematic studies on the effects 
of inclination and finite temperatures will be published 
elsewhere | p7| . 

As already mentioned, our results have important im- 
plications for growth processes on copper surfaces vici- 
nal to the (100) plane. First, under molecular beam epi- 
taxy (MBE) conditions step crossing should be very rare. 
Thus, our data can be used to construct a growth model 
for individual steps. In the case of copper, such a growth 
model ijl^ should include the following features: adatom 
motion along a straight step edge which is very fast, and 
motion through a kink site at the step edge which in 
turn is much slower than simple diffusion on a flat ter- 
race. In the simplest approximation, step crossing and 
evaporation can be neglected. The activation energies for 
these processes are well separated from each other and 
well below those of the neglected ones, which makes the 
model simple and hopefully applicable to a variety of vic- 
inal fee metal surfaces. This is supported by experiments 
[^,^jlj] . For some fee metals step crossings could be more 
significant, and then the detailed interplay between dif- 
ferent mechanisms at step edges and kink statistics has 
to be taken into account. Note that during growth a 
non-vanishing kink concentration is naturally provided 
by adatoms deposited on the terraces. The growth rules 
suggested above differ substantially from those expected 
to describe MBE on stepped surfaces of silicon , where 
the strong anisotropy of diffusion pi] results in an inter- 
esting step morphology ||5|,^. 
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FIG. 1. Ideal fcc(119) surface: (a) Perspective view and 
(b) top view. The size of the unit cell is shown by a dashed 
line. 



FIG. 2. The adiabatic potential experienced by an adatom 
on the Cu(119) surface as given by EMT at zero temperature: 
(a) a contour plot of the potential, where the global minimum 
is at (0,0) and the energy difference between each contour is 
0.1 eV, and (b) minimum energy route of an adatom diffusing 
in the x direction. 

FIG. 3. Snapshots of two exotic diffusion events from MD 
simulations of the Cu(119) surface (a) at T = 700 K and (b) 
at T = 750 K. These events only include surface atoms which 
are colored black. Atoms in the adjacent layers are shown in 
grey. Only part of the simulation cell is shown. In (a), the 
surface atom marked with a cross is pushed up and a hole is 
left behind at the descending step edge. In (b), the surface 
atom marked with a cross is pushed up and the vacancy left 
behind diffuses to the lower step on the right. 

FIG. 4. Dominant diffusion mechanisms at the step edge 
on an fcc(llm) surface (top view). Black circles are adatoms 
and open circles denote surface atoms. A shaded circle shows 
the position of an atom after the diffusion event has taken 
place. Kl is the mechanism with the lowest activation bar- 
rier of the several possibilities for hopping across a step near a 
kink. The five stages of the process X3 are shown by numbers 
in parentheses. The processes seen in Fig. 3(a) and 3(b) cor- 
respond to X2 and X3, respectively, albeit without an initial 
adatom on the lower terrace. 

TABLE 1. Activation energies of some diffusion mecha- 
nisms near step edges on Cu(llm) surfaces as given by EMT. 
The whole system around the adatom was allowed to relax 
in the calculations. The labels in the first column are those 
used in Fig. 4. The third column shows the barriers of the 
corresponding reversed processes. For the mechanism X3 the 
height of the highest barrier (stage 1) is given. 



M (mechanism) 


Em (eV) 


E^M (eV) 


HI 


0.867 


0.578 


H2 


0.258 


0.258 


XI 


0.909 


0.631 


X2 


1.310 


1.310 


X3 


1.074 




Kl 


0.842 


0.573 


K2 


0.782 


0.492 


K3 


0.879 


0.143 


K4 


0.518 


0.239 
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